Optimal quantum pump in the presence of a superconducting lead 
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We investigate the parametric pumping of a hybrid struc- 
ture consisting of a normal quantum dot, a normal lead and 
a superconducting lead. Using the time dependent scattering 
matrix theory, we have derived a general expression for the 
pumped electric current and heat current. We have also de- 
rived the relationship among the instantaneous pumped heat 
current, electric current, and shot noise. This gives a lower 
bound for the pumped heat current in the hybrid system sim- 
ilar to that of the normal case obtained by Avron et al. 
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Since n the seminar 

work of Brouwerjil the physics of|-paxametric pumping 
has attracted increasing attentionHiiJ Recently, Avron 
et allia have considered the pumped heat current in the 
adiabatic regime and found a general lower bound for 
the heat current. This defines an optimal pump if the 
heat current equals to the ppwer of Joule heat dissipated 
during the pumping processed. As a consequence, the op- 
timal pump is noiseless and charge transported is quan- 
tized. The physics of pumped heat curnpjt has also been 
investigated by Moskalets and Buttikerlia who have de- 
rived a general formula for the heat current in the weak 
pumping regime and the shot noise generated during the 
pumping process. In the strong pumping regime, the 
heat current has been studied within the time-dependent 
scattering matrix theoryEj and the existence of optimal 
pump has beaa. examined. For chaotic quantum dots, 
Polianski et alEll have developed a time-dependent scat- 
tering matrix theory to account for the shot noise for 
parametric pumping and mesoscopic fluctuation for arbi- 
trary temperature and beyond bilinear response. In this 
paper, we investigate the pumped heat current for a nor- 
mal superconducting (NS) hybrid system which consists 
of a normal quantum dot, a normal lead and a super- 
conducting lead. In the adiabatic regime, the energy of 
charged carriers (electron or hole) is within the super- 
conducting, energy gap and hence physics of the Andreev 
reflcctionEa dominates. We have derived a general expres- 
sion for the pumped electric current and heat current in 
the presence of superconducting lead which is valid at 
finite pumping amplitude and finite temperature. Our 
theorypis based on the time-dependent scattering matrix 
theoryB. Since our theory is perturbative in nature, go- 
ing beyond the adiabatic regime, one can in principle 
obtain the pumped electric current and heat current to 
any order in frequency. In the adiabatic regime, we have 
also derived a relationship among the instantaneous heat 



current, electric current, and the shot noise. This sets 
a lower bound forJdie pumped heat current. Similar to 
the normal systemEfl, a quantum pump will be optimal if 
the pumped heat current reaches its lower bound. As a 
result, the charge transported will be quantized and the 
system is noiseless just like the normal system. We have 
also compared with the pumped heat current of NS struc- 
ture with that of normal structure. For a single pumping 
potential, the total heat currents generated are the same 
for NS and normal systems. For two pumping potentials 
the pumped heat current for NS system can be either 
larger or smaller than that of normal system depending 
on the phase difference between two pumping potentials. 

For the purpose of presentation, we consider the 
pumped electric current first. We start with the general 
definition for the electric current of type a (electron or 
hole) in the .left normal lead in scattering matrix theory 
(ft = 5 = 1)^ 
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where < ... > denotes the quantum average and I e ,La is 
the electric current operator of type a in the left lead, 



g a [b{ a (t)b La (t) - a[ a (t)a La (t)} 



(2) 



Here the operators 6l q and a^a are annihilation oper- 
ators for the outgoing and incoming carriers of type a 
in the left lead and q a = 1,-1 i fof p = e i^- They are 
related by the scattering matrixlljlri] 



b La (t) =J2j dt'S af }(t,t')a L /s(t') 



(3) 



where the timc-dcpcndence of the scattering matrix S is 
due to the slowly time- varying pumping potential X(t). 
The distribution function can be obtained by taking the 
quantum aver age, lj 



< al a {E)a Lp (E') >= 5 af3 S(E - E')f L (E) 
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where ar, a (E) is the Fourier transform of a.La{t) and 
Jl{E) is the Fermi distribution function of the left lead. 
From Eqs.(||), (^), and (^), the pumped electric current 
is given by 

Qa f At f 
I e .La= lim / dt dt 1 dt 2 'y]S a p(t > t 1 ) 
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where f(t) = J(dE/2ir)exp(—iEt)f(E). After changing 
of the variable to = (t\ + ta)/2 and r = t\ — ti and 
using the following Wigner transform for the scattering 
matrixEl, 

S{t,t')=J d ^e-^-^S{E^) (6) 
then Eq.(|) becomes, 
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Changing the variables again to n = t — t and i' = 
(t + io)/2 and integrating over ri, we obtain 



Ie,La — 



lim [ dt'drdE e tET f(T) 
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J2S a p(E,t' + T/4)S* ap (E,t>-T/4)- qa J (8) 





Using the fact that 
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Eq.(H) becomes 
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where T p is the period of the pumping cycle and S is 
a 2 x 2 scattering matrix for NS structure with matrix 
element S a p (a,j3 = e,h). Eq.([lO|) is symbolic and is the 
central result of this paper. One can in principle obtain 
the pumped electric current to any order in frequency. 
For instance, to get the electric current up to oj, it is 
enough to expand f(E + id t /2) up to the first order in 
dt, from which we obtain 
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Note that from the unitary condition of the scattering 
matrix S, we have 



^ $af3 Sa{) ~ 1 





(12) 



taking the derivative with respect to time, we obtain 

~) d t S* a/3 S af3 + c.c. = (13) 





Hence Im[(9 t 5 t 5) cea ] i 
batic regime, we have 
dxSapdtXi + ...] where X = dX/dt. Up to the order 
u>, we can neglect the contribution from d x .s a p. At zero 
temperature, Eq.(0) becomes, 



— i(dtS^S) aa . In the adia- 
dtS af 3 — Yl l i[dx i S a pdtX l + 
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which aspees with the theory of nonequilibrium Green's 
functioned. 

Now we proceed to derive the pumped heat current for 
NS structure. We note that the heat current is defined 
as the particle current multiplied by the energy measured 
from the Fermi level, we thus have from Eq.dTT 
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Expanding the heat current up <9 t 2 and after some algebra, 



we finally obtained the heat current up to u> , 
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dt / dE d E f{d t S d t &) aa (16) 



Now we derive the relationship between instantaneous 
electric current (denoted as I e (t) = —iq a (dtS' i S) aa ) and 
heat current (I q (t) = (dtSWtS) aa ). Now the instanta- 
neous heat current becomes 



/,(*)= (d t & d t S) aa = (d t & dtS) c 







(17) 



We see that the diagonal term in Eq.(|17|) is just the I^(t) 
and the off diagonal term is t he . , s h ot noise So (t) generated 
during the pumping processt3£3E3. Therefore, we have 
the following relationship 

I q {t) = Il(t) + S (t) (18) 

or general lower bound for the heat current 

is) > m (19) 

The condition of optimal pump for NS .structure is de- 
fined as So{t) = 0. Following Avron et alia it is straight- 
forward to show that the charge transported through the 
system per cycle is quantized if the quantum pump is 
optimal. 

Now we consider the weak pumping limit for a sym- 
metric double barrier structure in the presence of super- 
conducting lead. The double barrier structure is modeled 
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by potential V(x) — X 1 (t)S(x + a) + X2(t)S(x — a) where 
Xi(t) = X a +X 1 sin(cji) andA 2 (t) = X +X 2 sin(u;i+0). 
In this limit, it is easy to show that the pumped electric 
current and heat current are given, respectively, by 



ujq a sin 4>X 1 X 2 
2^ 



Im[{d Xl &d x J) aa \ (20) 



and 



i q ,La = ^[xld Xl &d Xl s + x 2 2 d X2 &d x J 

+ 2cos(f>XiX 2 Re{d Xl Sd x J^] aa (21) 



In Eqs.(|2C|) and (|21j), we have set X, = in S after the 
partial derivatives. Now we will calculate the pumped 
heat current I q ^La=e for the double barrier NS system. 
For the NS rSje^em, the scattering matrix S ee and She 
are given by 



S — Su + (§12(1 — R1S22) Ri§2i 



where 



Sij(E) = 



Sij(E) 







(22) 



(23) 



with Sij being usual scattering matrix for the normal 
structure. i?j = aa x is the 2x2 scattering matrix at 
NS interface with off diagonal matrix element a. Here 
a = {E - WA 2 - E 2 )/A with v = 1 when E > -A 
and v = — 1 when E < — A. In Eq. (|22]) , the energy 
E is measured relative to the chemical potential fi of 
the superconducting lead. Eq.(22) has clear physical 
meaning 



The first term is the direct reflection from 
the normal scattering structure and the second term can 
be expanded as S12R1S21 + S12R1S22R1S21 + ■■■ which 
is clearly the multiple Andreev reflection in the hybrid 
structure. From Eq.(^2j) we obtain the well kaown ex- 
pressions for the scattering matrix S ee and ShJ^ 

S ee (E) = Sn(E) + a 2 Si2{E)S*2 2 {-E)M e S2i{E) (24) 

and 



S he (E) = aS{ 2 {-E)M e S 2 i{E) 



(25) 



with M e = [1 - a 2 S 22 (E)S^ 2 {-E)]- 1 . In the case of 
parametric pumping, we assume that the Fermi energy 
is in line with the chemical potential of superconducting 
lead, so E = and a — —i. For the symmetric NS system 
at resonance, we have Sn = and Si 2 = e~ 2tka in the 
absence of pumping potential. Therefore, from Eqs.(|2T 
and (|25|), we have 



and 



d Xl/2 S ee — d x Su - S 2 2 d x ,,,5*! 



d x S he = -i(d x Sl 2 Si2 + cc.) 
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(27) 



where we have used theJact that d Xl S2 2 = dx 2 Sn. Us- 
ing Fisher-Lee relationEZI S a p — + wG^ and the 
Dyson equation d Xj G r af} = G r aj G r j/3 c3 : we have d Xl Su = 
ivGi X G\i — —i/v, d X2 Su — ivG^G^i = —iS 2 2 /v, and 



d Xl „Su 



-1S12/V with the velocity v = 2k. Thus 



from Eqs.(|26j) and (|27 

d X She = 0. From Eq 



we have d Xj S e 
2l|), we obtain 



2d Xj Sn and 
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[X 2 + X 2 + 2 cos (f> X X X 2 cos 4fca] (28) 



which should be compared with the pumped heat current 
in the normal case 



%l = fin = £jp[ X i + X l + cos ^ X i^(l + cos4fca)] 

(29) 

We note that in the NS system, the heat current flows 
out only through the normal lead while for normal sys- 
tem, the heat current pumps out through both leads. 
Comparing Eqs.(|2^) and (p9|), we have 
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Hence the total pumped heat current generated in the 
normal system is can be either larger or smaller than 
that in the NS system depending on the phase difference 
of two pumping potentials. For a single pump, by set- 
ting X2 — in Eqs.(p8|) and (p9|), we see that the total 
pumped heat currents are the same for both NS and nor- 
mal systems. This is different from the pumped electric 
current where in the weak pumping regime at resonance, 
the electric current for-NS system is four time larger than 
that of normal systeml. 

In summary, we have derived a general expression for 
the pumped electric current and heat current in the pres- 
ence of superconducting lead using the time-dependent 
scattering matrix theory. Our theory is valid at finite 
pumping amplitude and can be applied to the multi- 
channel systems. Since our theory is perturbative in na- 
ture, we can expanding Eq.([T^) to the higher order in fre- 
quency and hence approach to the nonadiabatic regime. 
Our theory can also be easily extended to the case of 
multi-terminal structures. Although our expression is 
derived for NS system, it is also valid the normal sys- 
tem as well by simply replacing the NS scattering matrix 
S a p with a[3 — e, h by normal system scattering matrix 
Sij with i,j = 1,2 in Eqs.© and fiufy. For the NS 
system, we have found the lower bound fox.the pumped 
heat current similar to that of Avron et alE£l for the nor- 
mal system. As a result, the optimal pump can exist for 
NS system as well. In the weak pumping limit, we have 
examined the pumped heat current for NS system for a 
double barrier structure at resonance. For two parameter 
pump, we found that the total pumped heat current for 
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NS structure can be larger or smaller than that of normal 
structure depending on the phase different between two 
pumping parameters. 
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